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Abstract 



We present a method to measure the weak phase 7 in the three-body decay of 
charged mesons to the final states DK^tt . These decays are mediated by 
interfering amplitudes which are color-allowed and hence relatively large. As a 
result, large CP violation effects that could be observed with high statistical 
significance are possible. In addition, the three-body decay helps resolve 
discrete ambiguities that are usually present in measurements of the weak 
phase. The experimental implications of conducting these measurements with 
three-body decays are discussed, and the sensitivity of the method is evaluated 
using a simulation. 



1 



Typeset using REVTpjX 



I. INTRODUCTION 



CP violation is currently the focus of a great deal of attention. Since the start of oper- 
ation of the B-factories, the standard model description of CP violation via the Cabibbo- 
Kobayashi-Maskawa (CKM) matrix [1] is being tested with increasing precision. BaBar [2] 
and Belle [3] have recently published measurements of the CKM parameter sin(2/3), where 
(3 = arg (— VcdV'b/VtdVfo), verifying the CKM mechanism to within the experimental sen- 
sitivity. Although improved measurements of sin (2(3) and B s — B s mixing will probe the 
theory with greater scrutiny during the next few years, the measurement of the other angles 
of the unitarity triangle are necessary for a comprehensive study of CP violation. 

Important constraints on the theory will be obtained from measurements of the CKM 
phase 7 = arg(— V u dV* b /V c dV* b ). A promising method for measuring this phase in the B 
system has been proposed [4]. Although this method involves color- allowed decays and 
hence offers favorable rates, it makes use of B s mesons, which are not produced at B-Factories 
operating at the T(4S) resonance. By contrast, the extraction of 7 using the B M and B d 
system generally involves decays which are highly suppressed or difficult to reconstruct. In 
addition, these methods are generally subject to an eight-fold ambiguity due to a-priori 
unknown strong phases [4,5]. As a result, obtaining satisfactory sensitivity requires very 
high statistics and necessitates the use of as many decay modes and measurement methods 
as possible. 

One important class of theoretically clean measurements will make use of decays of the 
type B — > DK. Gronau and Wyler [6] have proposed to measure sin 2 7 in the interference 
between the b — > cus decay B + — > D°K + and the color-suppressed b — > ucs decay B + — > 
D°K + . Interference between these amplitudes takes place when the D meson is observed as 
one of the CP-eigenstates 



which are identified by their decay products, such as K + K~ or K s tt°. Several variations of 
this method have been developed [7-10], including addressing the effects of doubly Cabibbo- 
suppressed decays of the D meson [11] and mixing and CP-violation in the neutral D meson 
system [12], as well as insights to be gained from charm factory measurements [13]. 

A serious difficulty with measuring 7 using B — > DK is that the b — > ucs amplitude 
is expected to be extremely small. To a large degree, this is due to the color-suppression 
associated with the internal spectator diagram through which this amplitude proceeds. In 
the factorization model, color-suppression of the amplitude is parameterized by the phe- 
nomenological ratio |a 2 /ai|. This ratio is measured to be about 0.25 [14] by comparing 
decay modes which depend only on color-allowed amplitudes with those that depend on 
both color- allowed and color-suppressed amplitudes. With this value of 1 02/01!, one expects 
the amplitude ratio \A(B + -> D°K + )/A(B + -> WK + )\ to be only about 0.1. The small 
branching fraction B(B + — > D°K + ) is therefore very difficult to measure with adequate 
precision, resulting in a large statistical error in the measurement of sin 2 7. The recent 
observation of the color-suppressed decays B° — > ZX*)°7r°, D°i], and D°uj by Belle [15] and 
CLEO [16] has raised the possibility that (02/01 1 may be effectively larger in some modes. 
However, significant suppression is still expected for internal spectator diagrams. 
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This difficulty has led to attempts to address the problems presented by color suppression. 
Dunietz [7] proposed to apply the method to the decays B° — > DK*°, making use of the 
fact that the decay K*° — > K + ir~ tags the flavor of the B°. In this mode, both the b — > cos 
and the b ucs amplitudes are color suppressed, and hence of similar magnitudes, albeit 
small. Jang and Ko [9] and Gronau and Rosner [10] have devised a method in which the 
small branching fraction of the color-suppressed decay B + — > D°K + does not have to be 
measured directly. Rather, it is essentially inferred by using the larger branching fractions 
of the decays B° -> D-R+, B° -> WK° and B° -> D h2 K°. Quantitative study suggests 
that the various alternative methods are roughly as sensitive as the method of Ref . [6] , and 
are thus useful for increasing statistics and providing consistency checks [5]. 



II. MEASURING 7 WITH COLOR- ALLOWED B -> DKir DECAYS 

In this paper we investigate a way to circumvent the color suppression penalty by using 
B ± decay modes which could potentially offer significantly large branching fractions, as well 
as large CP asymmetries. Similar modes involving neutral B decays can also be used. For 
example the final state D (D )K ± tt t can be analyzed with the same technique as described 
here. Some other decays (such as B° — > D~K s tt + ) need a different treatment and will be 
discussed elsewhere [17]. The particular decays which are considered here are of the type 
B — > D^K^n(p). These three body final states may be obtained by popping a qq pair in 
color allowed decays. Although modes where one or more of the three final state particles 
is a vector can also be used, for clarity and simplicity only the mode B ± — > D ^^ is 
discussed here. 
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FIG. 1. Feynmann Diagrams for the de- 
cay B + -> WK + tt q involving the CKM ma- 
trix element product V* b V us 
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FIG. 2. Feynmann Diagrams for the de- 
cay B + -> D°K + ir° involving the CKM ma- 
trix element product V* b V cs 



Figs. 1 and 2 show the diagrams leading to the final states of interest. As can be seen, 
the leading diagrams (Fig. la and 2a) are both color-allowed and of order A 3 = sin(# c ) in the 
Wolfenstein parameterization [18], where 8 C is the Cabibbo mixing angle. Due to the absence 
of color suppression, both interfering amplitudes are large, avoiding the complications which 
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arise due to the small magnitude of the b —>■ ucs amplitude in the two-body decays. As a 
result, observable CP- violating effects in the three-body decays are expected to be large, and 
the b — > ucs amplitude is more easily measured from the relatively large branching fraction 
B{B + — > D°K + 7r°), which is now subject to significantly less contamination from doubly 
Cabibbo-suppressed D meson decays than the corresponding two-body modes. However, 
should the b — > ucs amplitude be unexpectedly small, one could still carry out the analysis 
described in this paper by taking doubly Cabibbo-suppressed decays into account [5,11]. 
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FIG. 4. Illustration of the triangle re- 
lations in the decays — > D^K^n and 
B ± -> D K±ir . 



nr(DV) 

FIG. 3. Two points on the Dalitz 
plot of the decays B + -> D^K+tt and 
B + _» D°K + ir°. 



Let us examine how one could observe CP violation and measure the angle 7 in these 
decays. We first consider the case of very large statistics, and then discuss how one would 
proceed when the data sample is limited. Since we are dealing with a three-body decay, we 
use the Dalitz plot of the system DK ± n° (see Fig. 3). Selecting a particular point % in this 
representation, Eq. (1) implies the relations 
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Let us write the amplitudes corresponding to the transitions in Fig. 1 and 2 as 



Ai(B + -> D°K + n°) = A Ci e i&c < , Ai{B+ -> D°K^ 
Ai(B~ -> D°K-n°) = A Ci e iSc > , Ai(B~ -> WK~ 
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(3) 



where 7 is the relative phase of the CKM matrix elements involved in this decay, and Ac 
{Ay) and 5c (Su) are the real amplitude and CP-conserving strong interaction phase of 
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the transitions of Fig. 1 (Fig. 2). Let us note here that we have used the Wolfenstein 
parameterization at order 0(X 3 ). Should one use the full expansion, a small weak phase of 
order A 4 would be present. Indeed the angle which is measured is 7' = aig(V u dV* b /V cs V* b ) = 
7 + £, where £ = arg (— V cd V* s /V ud V* s ). The angle £ is one of the angles [19] arising from the 
unitarity relation V ud V* s + V c dV* s + V t dV t * s = 0. The amplitudes in Eqs. 3 can be obtained 
from the measurements of the B decay widths 

_> WK + ti q ) = Fi(B- -> D°K-n°) = A c \ 
F t (B + -> D°K + 7i°) = Fi(B- -> WK'tt ) = A u l (4) 

Eq. (2) implies 

2F Z (B + - Dl 2 K + 7T°) = Ad + A, 2 ± 2Ar*A^ cos(A^ + 7) 

2F i (B~ - D1 2 K-tt°) = Ad + ^2 ± 2^4* cos(A5, - 7), (5) 

where A5i = — 5c«- Thus, by measuring the widths in Eq. 4 and 5, one extracts sin 2 7 
from 

sin 2 7 = \ (l - C^±V(1-C 2 )(1-^)) , (6) 

where C = cos(A5j + 7) and C = cos(A5j — 7). Hence in the limit of very high statistics, 
one would extract sin 2 7 for each point i of the Dalitz plot and therefore obtain many 
measurements of the same quantity. This would allow one to obtain a large set of redundant 
measurements from which a precise and consistent value of sin 2 7 could be extracted. We 
note that in our treatment we disregard doubly Cabibbo-suppressed D° decays [11] since, 
due to the lack of color-suppression, their effect is small, and can be dealt with [5] in any 
case. 

In every point of the Dalitz plot, 7 is obtained with an eight-fold ambiguity, which is a 
consequence of the invariance of the cos(A<5j ±7) terms in Eq. (5) under the three symmetry 
operations [5] 

S cx : 7 -> 5 ,5^7 

S sig n : 7 -»■ -7 , $ -»■ S ( 7 ) 

Sit '■ 1 — >7 + 7r,5-^5 + 7T. 

However, an important benefit is gained from the multiple measurements made in different 
points of the Dalitz plot. When results from the different points are combined, some of the 
ambiguity will be resolved, in the likely case that the strong phase A5i varies from one region 
of the Dalitz to the other. This variation can either be due to the presence of resonances 
or because of a varying phase in the non-resonant contribution. In this case, the exchange 
symmetry S cx is numerically different from one point to the other, which in effect breaks 
this symmetry and resolves the ambiguity. 

Similarly, the .S^gn symmetry is broken if there exists some a-priori knowledge of the de- 
pendence of A5{ on the Dalitz plot parameters. This knowledge is provided by the existence 
of broad resonances, whose Breit-Wigner phase variation is known and may be assumed to 
dominate the phase variation over the width of the resonance. To illustrate this, let i and j 
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be two points in the Dalitz plot, corresponding to different values of the invariant mass of 
the decay products of a particular resonance. For simplicity we consider only one resonance. 
One then measures cos(A<5j ± 7) at point i and cos(A5j + ± 7) at point j, where 
is known from the parameters of the resonance. It is important to note that the sign of 
ctij is also known, hence it does not change under S S i gn . Therefore, should one choose the 
*Ssign-related solution cos(— A5i =F 7) at point i, one would get cos(— A<5« + =F 7) at point 
j. Since this is different from cos(A5j + ± 7), the S^gn ambiguity is resolved. This is 
illustrated graphically in Eq. 8: 

cos(A5j ± 7) cos(— A5i =f 7) 

BW I I BW (8) 

^sign 

cos(A5j + atij ± 7) ^-f-^ cos(-A5j + =p 7) 

Thus, broad resonances reduce the initial eight-fold ambiguity to the two-fold ambiguity 
of the S n symmetry, which is not broken. Fortunately, S n leads to the well-separated so- 
lutions 7 and 7 + 7r, the correct one of which is easily identified when this measurement is 
combined with other measurements of the unitarity triangle. 



III. THE FINITE STATISTICS CASE 

Since experimental data sets will be finite, extracting 7 will require making use of a 
limited set of parameters to describe the variation of amplitudes and strong phases over 
the Dalitz plot. The consistency of this approach can be verified by comparing the results 
obtained from fits of the data in a few different regions of the Dalitz plot, and the systematic 
error due to the choice of the parameterization of the data may be obtained by using different 
parameterizations. 

A fairly general parameterization assumes the existence of Nr Breit-Wigner resonances, 
as well as a non-resonant contribution: 

Az(B + -> WK + n°) = (ac e i5c ° + J^A Cj B s ^) e^j e iSc ® 

At(B + - D°K + n°) = (a U0 e i5 "° + A U] B Sj (0 e^j e^e l \ (9) 
where £ represents the Dalitz plot variables, 

II SO l>, iO<'* - (10) 

is the Breit-Wigner amplitude for a particle of spin Sj, normalized such that / (b s (£)) 2 <i£ = 1, 
Auq and 5u {Acq and Sco) are the magnitude and CP-conserving phase of the non-resonant 
b — > ucs (b — > cus) amplitude, and Auj and 5uj (A C j and 5c j) are the magnitudes and 
CP-conserving phase of the b — * ucs (b — > cus) amplitude associated with resonance j [20]. 
The functions 8c (£) and 5u(£) may be assumed to vary slowly over the Dalitz plot, allowing 
their description in terms of a small number of parameters. Eq. (1) again implies 
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A^(B + -> Dl 2 K+n°) = -L(^(£?+ -> D°K + n°) ± ^(5+ -> ZW+vr )). (11) 

The decay amplitudes of B _ mesons are identical to those of Eqs. (9) and (11), with 7 
replaced by —7. 

The decay amplitudes of Eqs. (9) and (11) can be used to conduct the full data analysis. 
This is done by constructing the probability density function (PDF) 

p(0 = MW)I 2 , (12) 

where the amplitude A%(f) is given by one of the expressions of Eq. (9), Eq. (11), or their 
CP-conjugates, depending on the final state /. Given a sample of N e signal events, 7 and 
the other unknown parameters of Eq. (9) are determined by minimizing the negative log 
likelihood function 

X 2 = -2$>gP(&), (13) 
i=\ 

where & are the Dalitz plot variables of event i. 



IV. RESONANCES AND AMBIGUITIES 

It is worthwhile to consider the resonances which may contribute to the D°ir ± 7r° final 
state. Obvious candidates are broad D** and D** states. However, only the ones which 
can decay as D**° — > D°n° or D** + — > D°K + are relevant for the final state of interest. 
This exclude the 1 + states, which would decay to D*n or D*K. Furthermore, since the 
D** + is essentially produced through a W + , the 2 + state is forbidden as well. Thus, one 
does not expect a large contribution from these states. A promising candidate might be the 
broad Dq° recently observed by the Belle collaboration [21]. We note that including such 
resonances in the analysis does not raise particular difficulties and would further enhance 
the sensitivity of the 7 measurement. Similar arguments can be made for higher excited K 
states. 

One also expects narrow resonances, such as the D*(2007)° and a narrow D** + state, to 
be produced. However, as seen in the Dalitz plot of Fig. 5, these resonances do not overlap, 
and hence do not interfere. In addition, interference between a very narrow resonance and 
either a broad resonance or a non-resonant term is suppressed in proportion to the square 
root of the narrow resonance width. Therefore, narrow resonances contribute significantly 
to the CP violation measurement only if both the b — > cus and b — > ucs amplitudes proceed 
through the same resonance. This scenario is favorable, but is not necessary for the success 
of our method, and will therefore not be focused on in the rest of this study. 

In what follows, we discuss important properties of the method by considering the illus- 
trative case, in which the b — > ucs decay proceeds only via a non-resonant amplitude, and 
the b — > cus decay has a non-resonant contribution and a single resonant amplitude. For 
concreteness, the resonance is taken to be the K* ± (892). We take the ^-dependent non- 
resonant phases to be 5c(£) = 8u(0 — 0. Under these circumstances, the PDF of Eq. (12) 
depends on four cosine terms that are measured in the experiment: 
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c oo 

„± 

-K*0 



cos(5 m ± 7) 

cos(5 W0 - S K > - 6 K *(g) ± 7), 



(14) 



where 5 K *(£) is the ^-dependent K* Breit-Wigner phase of Eq. (10). The cosines Cq (c^* ) 
arise from interference between the non-resonant (resonant) b — > cus amplitude and the 
non-resonant b — > ucs amplitude. 

The phases 5uo, 5 k*, and 7 are all a-priori unknown. However, it is important to note 
that 8k* is fully determined from the interference between the resonant and non-resonant 
contributions to the relatively high statistics decay mode B + — > D K + n° as a function of 
the Dalitz plot variables. Therefore, 5k* is obtained with no ambiguities, and with an error 
much smaller than those of 5uo or 7*. Consequently, the only relevant symmetry operations 
are 



'sign 



oK* + 
^cx 

CK*- 



7 

7 
7 
7 
7 



$uo 
-7 

7 + 7T 

duo ~ $K* 
—5uo + 5k* 



$U0 

5uo 
5uo 

$U0 



7 

-5 m 

5uo + 7T 
7 + 8k* 

— 7 + $K* 



(15) 



As discussed above, only S n is a symmetry of all four cosines of Eq. (14), and is therefore 
fully unresolved. The transformation properties of the cosines under any combination of the 
remaining four operations that can lead to an ambiguity are shown in Table I. 



TABLE I. Invariance of each of the cosines of Eq. (14) under combinations of the symmetry 
operations of Eq. (15), excluding S n . Full invariance (approximate invariance) is indicated by a 

vV)- 



Operation 


4*o 


C K*0 


°00 


c 00 






Non-resonant regime 






Sex 


V 




V 


V 


5sign 






V 


V 


"Sex S'sign 




V 


V 


V 






Resonant regime 






qK* + 


V 


7 


V 




oK*- 
D ex 


V 


V 




V 


qK*+ qK* — 
°ex °cx 


V 


V 







*Even when one of the b — > cus amplitudes is small enough that the determination of 5k* becomes 
difficult, one effectively has 5uo —> 5uo ~ 5k*, 5k* — ► 0, and the determination of 5k* is again not 
a problem. 



8 



While none of the operations leaves all four cosines invariant, it is important to note 
cases where c K * are approximately invariant under S^* + , S^.*~, or their product. We 
define approximate invariance under the operation S to be 

Sap P c|*o(<MO) = c|*o(-<MO)- ( 16 ) 

Approximate invariance arises due to the fact that far from the peak of the K* resonance, 
5k* (0 changes slowly as a function of the Kn invariant mass, and takes values around and 
7r. Therefore, for events in the tails of the Breit-Wigner, 5k* (0 is almost invariant under 
any S'app satisfying Eq. (16). One can see that approximate invariance of one of the cosines 
c K * implies minimal change in the x 2 of Eq. (13), which may result in a resolved yet clearly 
observable ambiguity. Since both c K *o terms are only approximately invariant under the 
product S^* + S^*~ , this ambiguity is more strongly resolved than either S^* + or S^*~ . 

Observing that no single operation in the Table I is a good symmetry of all cosines, 
one identifies two different regimes: In the non-resonant regime, interference with the non- 
resonant b — > cus is dominant, and only S cx and S S i gn may lead to ambiguities. In the 
resonant regime, the K* amplitude strongly dominates the b — > cus decay, and S^* + and 
S^*~ become the important ambiguities. In the transition between these regimes, the op- 
erations of Table I do not lead to clear ambiguities, as we have verified by simulation (See 
sec. V). Thus, while naively one may expect a 2 5 -fold ambiguity, in practice the observable 
ambiguity is no larger than eight-fold, with only the two-fold S n being fully unresolved, in 
the likely case of non-negligible resonant contribution. This is demonstrated in Fig. 8. Fur- 
thermore, although one may write down more products of the operations S cx , S sign , S^* + , 
and S^*~, only the products listed in Table I result in full or partial invariance of both 
cosines which dominate the same regime. The additional products do not result in any 
noticeable ambiguities. 



V. MEASUREMENT SENSITIVITY AND SIMULATION STUDIES 

To study the feasibility of the analysis using Eq. (13) and verify the predictions of Sec. IV, 
we conducted a simulation of the decays B ± — > D°K ± rr , — > D°K ± tt°, and B ± — ► 
Dia^Tr . Events were generated according to the PDF of Eq. (12), with the base parameter 
values given in Table II. In this table and throughout the rest of the paper, we use a tilde to 
denote the "true" parameter values used to generate events, while the corresponding plain 
symbols represent the "trial" parameters used to calculate the experimental \ 2 ■ 

The only non-vanishing amplitudes in the simulation were the non-resonant amplitudes 
in the b — > cus and b — > ucs decays, and the K* resonant b — > cus amplitude. For simplicity, 
additional resonances were not included in this demonstration. However, broad resonances 
that are observed in the data should be included in the actual data analysis. 

The simulations were conducted with a benchmark integrated luminosity of 400 fb _1 , 
which each of the asymmetric B-factories plan to collect by about 2005. The final state 
reconstruction efficiencies were calculated based on the capabilities of current T(4S) detec- 
tors. We assumed an efficiency of 70% for reconstructing the K ± , including track quality 
and particle identification requirements, and 60% for reconstructing the tt°. The product of 
reconstruction efficiencies and branching fractions of the D°, summed over the final states 
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TABLE II. Parameters used to generate events in the simulation. The value of Aqk* is chosen 
so as to roughly agree with the measurement of the corresponding branching fraction [22], taking 
into account the K* + — ► K + ir° branching fraction. 



Parameter 


Value 


Parameter 


Value 


7 


1.20 


Auo/Aco 


0.4 


~Sc(0 = s u (0 





A C K*/Aco 


1.0 


5k* 


1.8 


A-CK* 


~ V2 x io- 4 r B 




0.4 







TABLE III. The numbers of events obtained by averaging 100 simulations using the parameters 
of Table II and the reconstruction efficiencies listed in the text. 

Mode Signal events per 400fb~ 1 

B+ Wk+tt° = B~ -► D°K~7r° 2610 

B + D°K+7r° = B~ -► IfiK-Tr 205 

B + D h2 K + 7r° 186 

B~ -► Di 2^ _ vr° 234 



K~n + , K~7T + 7T°, and i^~7r + 7r~7r + , is taken to yield an effective efficiency of 6%. Using the 
CP-eigenstate final states K + K~, 7i + n^, K S 7f°, and K s p°, the effective efficiency for the 
sum of the D x and D 2 final states is 0.8%. All efficiencies are further reduced by a factor of 
1.7, in order to approximate the effect of background. The numbers of signal events obtained 
in each of the final states with the above efficiencies and the parameters of Table II are listed 
in Table III. 

In Figs. 6 through 8, we show the dependence of x 2 on the values of 7 and 5uo- The 
smallest value of x 2 is shown as zero. At each point in these figures, x 2 is calculated with the 
generated values of the amplitude ratios A Uo /Ac = A U q/Ac and A ck */A C q = A CK */Ac . 
We note that when these amplitude ratios are determined by a fit simultaneously with the 
phases, the correlations between the amplitudes and the phases are generally found to be 
less than 20%. Therefore, the results obtained with the amplitudes fixed to their true values 
are sufficiently realistic for the purpose of this demonstration. 

For each of these figures, we also show the one-dimensional minimum projection x 2 {l) = 
min{x 2 (7, Su )}, showing the smallest value of x 2 f° r each value of 7. 

Fig. 6 is a simulation obtained with the parameters of Table II, but with Aqk* — 0. 
With no resonant contribution, the eight-fold ambiguity of the perfect non-resonant regime 
is clearly visible. This would be the typical case for two-body final states. 

Fig. 7 is obtained with the parameters of Table II, but with A Co = 0. With no non- 
resonant b — > cus contribution, the eight-fold ambiguity of the perfect resonant regime is 
seen. The ambiguities corresponding to approximate invariance are clearly resolved, with 
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FIG. 5. Dalitz plots obtained from a simulation of B + and B~ decays into all final state, 
D ^^ , ZJOi^vr , and D 1:2 K ± 7r°. with the parameters of Table II. Along with non-resonant 
contributions, the resonances K*°, D*°, and D** are shown. 



the doubly- approximate ~ ambiguity resolved more strongly. 

Fig. 8 is obtained with the parameters of Table II and shows how efficient the method 
described in this paper could be for extracting the angle 7. With equal resonant and non- 
resonant b — > cus amplitudes, only the non-resonant regime ambiguities are observed, due to 
the relative suppression of the resonant interference terms discussed in Sec. IV. Nonetheless, 
the c^* terms are significant enough to resolve all but the S w ambiguity. S s - lgn is more 
strongly resolved, since it leaves neither of the c^* terms invariant. 

Also shown in Fig. 8 (dashed line) is the minimum projection of x 2 (l) °f the average 
experiment. This plot is obtained by averaging over x 2 {l) °f 50 simulated experiments, each 
generated with the parameters of Table II, but with different initial random numbers. The 
three standard deviation region of 7 allowed by the average experiment is indicated. This 
region spans the range [—0.06, 1.65], giving an idea of the sensitivity that may be obtained 
with these parameter values. 

In Figs. 9 through 12 we present <7 7 , the statistical error in the measurement of 7, obtained 
by fitting simulated event samples using the MINUIT package [23] , as a function of one of the 
parameters of Table II. All the other parameters were kept at the values listed in Table II. 
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FIG. 6. Top: x 2 as a function of 7 and Suo, with the parameters of Table II and no resonant 
contribution (Ack* = 0). Bottom: Minimum projection of x 2 onto 7. 



Each point in these plots is obtained by repeating the simulation 250 times, to minimize 
sample-to-sample statistical fluctuations. In all cases, all the parameters of Table II were 
determined by the fit. The arrows in these figures indicate the point corresponding to the 
parameters of Table II. The total number of signal events in all final states combined is the 
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FIG. 7. Top: x 2 as a function of 7 and <5^(b with no non-resonant b — ► cms contribution 
(^4co = 0)- The value <5j^* = 1.2 is used to ensure that ambiguities do not overlap. All other 
parameters are those of Table II. Bottom: Minimum projection of x 2 onto 7. 

same for each of the data points. The error bars describe the statistical error at each point, 
which is determined by the number of experiments simulated. 

One observes that <r 7 does not depend strongly on 8k*, and has a mild dependence on 
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FIG. 8. Top: x 2 as a function of 7 and <5^o> with the parameters of Table II. Bottom: The 
solid line shows the minimum projection of \ 2 onto 7 for the example experiment. The dashed line 
represents the average of 50 simulated experiments. The three standard deviation allowed range 
of 7 obtained from the average is indicated by arrows. 

5uo- As expected, strong dependence on Au / Acq is seen in Fig. 12. However, it should 
be noted that 07 changes very little for all values of Auq/Aqq above about 0.4, given that 
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FIG. 9. The error in 7, cr 7 , as a function of 7. 



the total number of signal events in all modes was kept constant in our simulation. This 
suggests that the likelihood for a significantly sensitive measurement is high over a broad 
range of parameters. With the parameters of Table II, we find a 1 ~ 0.23 = 13° with an 
integrated luminosity of 400 fb _1 . 



VI. CONCLUSIONS 

We have shown how 7 may be measured in the color-allowed decays B — > D^K^n(p), 
focusing on the simplest mode B^ — > D°K ± ir°. The absence of color suppression in the 
b — > ucs amplitudes is expected to result in relatively large rates and significant CP violation 
effects, and hence favorable experimental sensitivities. Although the Dalitz plot analysis 
required for this purpose constitutes some experimental complication, it should not pose 
a major difficulty, while being very effective at reducing the eight-fold ambiguities that 
constitute a serious limitation with other methods for measuring 7. Only the two-fold S n 
ambiguity cannot be resolved solely by our method, requiring additional constraints from 
other measurements of the unitarity triangle. As a result of these advantages, this method 
is likely to lead to relatively favorable errors and provide a significant measurement of 7, 
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FIG. 10. The error in 7, <r 7 , as a function of el- 



even with the current generation of B factory experiments. 
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